Cumulative Bar

The cumulative chart shows accumulation of both COVID-19 cases and Shooting cases in NYC county from 2020-03-01 until 2021-09-30. In order to make the plot more readable, We divided the COVID-19 cumulative cases by 1000 due to its rapid growth rate compared to shooting cases. The bar chart shows a potential peak shifting fluctuation, instead of co-frequency resonance, which may be due to the delayed effect of Covid-19 on society. Therefore, we go deeper into the relationship between the accumulation of shooting cases and COVID cases.

plot_ly(shooting_covid, x = ~date, y = ~covid_cum_per1000 , type = "bar", alpha = 1, name = "COVID cases(/1000)") %>% 
  add_trace(shooting_covid, x = ~date, y = ~shooting_cum, type = "bar", alpha = 0.6, name = "Shooting case") %>% 
  layout(
    title = "Cumulative bar for Covid and Shooting",
    barmode = "overlap",
    xaxis = list(title = "date"),
    yaxis = list(title = "cumulative scale")) %>% 
  layout(showlegend = TRUE, legend = list(font = list(size = 8)))

Linear relation between accumulations

The linear plot clearly shows a steep rise in shooting cases when COVID-19 case accumulated to 200,000 around 2020-05-15 in NYC. The growth rate becomes slower for nearly half a year from around 2020-11-15, and then increases again since 2021-05, but not as rapid as in 2020. Combined with the month and shooting cases analysis before, increase in shooting case during May is common, but especially high growth rate during the begining of epidemic is unusual. Whatโ€™s more, it is noteworthy that the spread of Delta virus also begins in the end of April, 2021.

shooting_covid %>% 
  mutate(text_label = str_c("date: ", date)) %>% 
  plot_ly(x = ~covid_cum_per1000, y = ~shooting_cum,type = 'scatter', mode = 'lines', text = ~text_label, alpha = 0.8, colors = "viridis") %>% 
  layout(
    title = "Linear Relationship between Accumulations ",
    xaxis = list(title = "Cumulative covid cases/1000"),
    yaxis = list(title = "Cumulative shooting cases"))
LS0tCnRpdGxlOiAiQ09WSUQiCm91dHB1dDoKICBodG1sX2RvY3VtZW50OgogICAgdG9jOiB0cnVlCiAgICB0b2NfZmxvYXQ6IHRydWUKICAgIGRmX3ByaW50OiBwYWdlZAogICAgY29kZV9mb2xkaW5nOiBoaWRlCiAgaHRtbF9ub3RlYm9vazoKICAgIHRvYzogeWVzCmVkaXRvcl9vcHRpb25zOgogIGNodW5rX291dHB1dF90eXBlOiBpbmxpbmUKLS0tCmBgYHtyIHNldHVwLCBpbmNsdWRlPUZBTFNFfQpsaWJyYXJ5KHRpZHl2ZXJzZSkKbGlicmFyeShwbG90bHkpCmxpYnJhcnkoZHBseXIpCmxpYnJhcnkodmlyaWRpcykKbGlicmFyeShnZ2FuaW1hdGUpCmxpYnJhcnkoZ2dtYXApCmxpYnJhcnkocG5nKQpsaWJyYXJ5KGdpZnNraSkKdGhlbWVfc2V0KHRoZW1lX21pbmltYWwoKSArIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbiA9ICJib3R0b20iKSkKCmtuaXRyOjpvcHRzX2NodW5rJHNldCgKICBmaWcud2lkdGggPSA4LAogIGZpZy5oZWlnaHQgPSA2LAogIG91dC53aWR0aCA9ICI5MCUiLAogIG1lc3NhZ2UgPSBGQUxTRSwKICB3YXJuaW5nID0gRkFMU0UsCiAgZXJyb3IgPSBGQUxTRQopCgpvcHRpb25zKAogIGdncGxvdDIuY29udGludW91cy5jb2xvciA9ICJ2aXJpZGlzIiwKICBnZ3Bsb3QyLmNvbnRpbnVvdXMuZmlsbCA9ICJ2aXJpZGlzIgopCgpzY2FsZV9jb2xvcl9kaXNjcmV0ZSA9IHNjYWxlX2NvbG9yX3ZpcmlkaXNfZApzY2FsZV9maWxsX2Rpc2NyZXRlID0gc2NhbGVfZmlsbF92aXJpZGlzX2QKCmBgYAoKYGBge3IgcmVhZCwgaW5jbHVkZSA9IEZBTFNFIH0KY292aWRfMTlfZGYgPSByZWFkLmNzdigiLi9kYXRhL2NvdmlkXzE5X255Yy5jc3YiKSAlPiUgCiAgc2VsZWN0KHllYXIsIG1vbnRoLCBkYXksIGNhc2VzKSAlPiUKICAgICBtdXRhdGUoCiAgICBtb250aCA9IGFzLmZhY3Rvcihtb250aCksIAogICAgZGF5ID0gYXMuZmFjdG9yKGRheSksIAogICAgeWVhciA9IGFzLmZhY3Rvcih5ZWFyKQogICApCgpgYGAKCgpgYGB7ciBpbXBvcnQgZGF0YSwgaW5jbHVkZSA9IEZBTFNFfQpoaXN0b3J5X2RmID0gcmVhZC5jc3YoJy4vZGF0YS9oaXN0b3JpYy5jc3YnKQoKeWVhcl90b19kYXRlX2RmID0gcmVhZC5jc3YoJy4vZGF0YS95ZWFyX3RvX2RhdGUuY3N2JykgJT4lIAogIHJlbmFtZShMb25fTGF0ID0gTmV3Lkdlb3JlZmVyZW5jZWQuQ29sdW1uKQpgYGAKCmBgYHtyIGRhdGEgY2xlYW5pbmcsIGluY2x1ZGUgPSBGQUxTRX0KU2hvb3RpbmdfZGYgPSBiaW5kX3Jvd3MoaGlzdG9yeV9kZiwgeWVhcl90b19kYXRlX2RmKSAlPiUgCiAgamFuaXRvcjo6Y2xlYW5fbmFtZXMoKSAlPiUgCiAgc2VwYXJhdGUob2NjdXJfZGF0ZSwgYygibW9udGgiLCJkYXkiLCJ5ZWFyIikpICU+JSAKICBzZXBhcmF0ZShvY2N1cl90aW1lLCBjKCJob3VyIiwibWludXRlIiwic2Vjb25kIikpICU+JSAKICBtdXRhdGUoCiAgICB5ZWFyID0gYXMubnVtZXJpYyh5ZWFyKSwgCiAgICBtb250aCA9IGFzLm51bWVyaWMobW9udGgpLAogICAgZGF5ID0gYXMubnVtZXJpYyhkYXkpLAogICAgaG91ciA9IGFzLm51bWVyaWMoaG91ciksCiAgICBtaW51dGUgPSBhcy5udW1lcmljKG1pbnV0ZSkKICAgICkgJT4lIAogIHNlbGVjdCgtc2Vjb25kKQpgYGAKCmBgYHtyIGRhdGEsIGluY2x1ZGU9RkFMU0V9CnNob290aW5nID0gU2hvb3RpbmdfZGYgJT4lIAogIHNlbGVjdCh5ZWFyLCBtb250aCwgZGF5KSAlPiUgCiAgZ3JvdXBfYnkoeWVhciwgbW9udGgsIGRheSkgJT4lIAogIHN1bW1hcmlzZShkYXljb3VudCA9IG4oKSkgJT4lIAogIHN1YnNldCh5ZWFyID09IDIwMjAgfCB5ZWFyID09IDIwMjEpCnN1YnNldDIgPSBzaG9vdGluZyAlPiUgCiAgc3Vic2V0KHllYXIgPT0gMjAyMCAmIG1vbnRoIDw9IDIpCnNob290aW5nID0gYW50aV9qb2luKHNob290aW5nLHN1YnNldDIpCgpzaG9vdGluZ19jb3ZpZCA9IG1lcmdlKGNvdmlkXzE5X2RmLCBzaG9vdGluZywgYnkgPSBjKCd5ZWFyJywgJ21vbnRoJywgJ2RheScpLCBhbGwgPSBUUlVFKSAlPiUgCiBtdXRhdGVfYWxsKH5yZXBsYWNlKC4sIGlzLm5hKC4pLCAwKSkgJT4lIAogbXV0YXRlKAogICB5ZWFyID0gYXMuY2hhcmFjdGVyKHllYXIpLCAKICAgbW9udGggPSBhcy5udW1lcmljKG1vbnRoKSwKICAgZGF5ID0gYXMubnVtZXJpYyhkYXkpLAogICBkYXRlID0gbHVicmlkYXRlOjptYWtlX2RhdGUoeWVhciwgbW9udGgsIGRheSkKICAgKSAlPiUgCiAgc2VsZWN0KGRhdGUsIGRheWNvdW50LCBjYXNlcykgJT4lIAogIHJlbmFtZSgKICAgIHNob290aW5nX2Nhc2VzID0gZGF5Y291bnQsCiAgICBjb3ZpZF9jdW0gPSBjYXNlcwogICkgJT4lIAogIG11dGF0ZSgKICAgIGRhdGUgPSBhcy5mYWN0b3IoZGF0ZSkKICApICU+JSAKICBtdXRhdGUoY292aWRfY2FzZXMgPSBjb3ZpZF9jdW0gLSBsYWcoY292aWRfY3VtLCBkZWZhdWx0ID0gMCwgb3JkZXJfYnkgPSBkYXRlKSkgJT4lIAogIG11dGF0ZSgKICAgIHNob290aW5nX2N1bSA9IGN1bXN1bShzaG9vdGluZ19jYXNlcyksCiAgICBjb3ZpZF9jdW1fcGVyMTAwMCA9IGNvdmlkX2N1bS8xMDAwCiAgICApCiAgCmBgYAoKIyMgQ3VtdWxhdGl2ZSBCYXIKClRoZSBjdW11bGF0aXZlIGNoYXJ0IHNob3dzIGFjY3VtdWxhdGlvbiBvZiBib3RoIENPVklELTE5IGNhc2VzIGFuZCBTaG9vdGluZyBjYXNlcyBpbiBOWUMgY291bnR5IGZyb20gMjAyMC0wMy0wMSB1bnRpbCAgMjAyMS0wOS0zMC4gSW4gb3JkZXIgdG8gbWFrZSB0aGUgcGxvdCBtb3JlIHJlYWRhYmxlLCBXZSBkaXZpZGVkIHRoZSBDT1ZJRC0xOSBjdW11bGF0aXZlIGNhc2VzIGJ5IDEwMDAgZHVlIHRvIGl0cyByYXBpZCBncm93dGggcmF0ZSBjb21wYXJlZCB0byBzaG9vdGluZyBjYXNlcy4gVGhlIGJhciBjaGFydCBzaG93cyBhIHBvdGVudGlhbCBwZWFrIHNoaWZ0aW5nIGZsdWN0dWF0aW9uLCBpbnN0ZWFkIG9mIGNvLWZyZXF1ZW5jeSByZXNvbmFuY2UsIHdoaWNoIG1heSBiZSBkdWUgdG8gdGhlIGRlbGF5ZWQgZWZmZWN0IG9mIENvdmlkLTE5IG9uIHNvY2lldHkuIFRoZXJlZm9yZSwgd2UgZ28gZGVlcGVyIGludG8gdGhlIHJlbGF0aW9uc2hpcCBiZXR3ZWVuIHRoZSBhY2N1bXVsYXRpb24gb2Ygc2hvb3RpbmcgY2FzZXMgYW5kIENPVklEIGNhc2VzLgpgYGB7ciBiYXJwbG90fQpwbG90X2x5KHNob290aW5nX2NvdmlkLCB4ID0gfmRhdGUsIHkgPSB+Y292aWRfY3VtX3BlcjEwMDAgLCB0eXBlID0gImJhciIsIGFscGhhID0gMSwgbmFtZSA9ICJDT1ZJRCBjYXNlcygvMTAwMCkiKSAlPiUgCiAgYWRkX3RyYWNlKHNob290aW5nX2NvdmlkLCB4ID0gfmRhdGUsIHkgPSB+c2hvb3RpbmdfY3VtLCB0eXBlID0gImJhciIsIGFscGhhID0gMC42LCBuYW1lID0gIlNob290aW5nIGNhc2UiKSAlPiUgCiAgbGF5b3V0KAogICAgdGl0bGUgPSAiQ3VtdWxhdGl2ZSBiYXIgZm9yIENvdmlkIGFuZCBTaG9vdGluZyIsCiAgICBiYXJtb2RlID0gIm92ZXJsYXAiLAogICAgeGF4aXMgPSBsaXN0KHRpdGxlID0gImRhdGUiKSwKICAgIHlheGlzID0gbGlzdCh0aXRsZSA9ICJjdW11bGF0aXZlIHNjYWxlIikpICU+JSAKICBsYXlvdXQoc2hvd2xlZ2VuZCA9IFRSVUUsIGxlZ2VuZCA9IGxpc3QoZm9udCA9IGxpc3Qoc2l6ZSA9IDgpKSkKYGBgCgojIyBMaW5lYXIgcmVsYXRpb24gYmV0d2VlbiBhY2N1bXVsYXRpb25zCgpUaGUgbGluZWFyIHBsb3QgY2xlYXJseSBzaG93cyBhIHN0ZWVwIHJpc2UgaW4gc2hvb3RpbmcgY2FzZXMgd2hlbiBDT1ZJRC0xOSBjYXNlIGFjY3VtdWxhdGVkIHRvIDIwMCwwMDAgYXJvdW5kIDIwMjAtMDUtMTUgaW4gTllDLiBUaGUgZ3Jvd3RoIHJhdGUgYmVjb21lcyBzbG93ZXIgZm9yIG5lYXJseSBoYWxmIGEgeWVhciBmcm9tIGFyb3VuZCAyMDIwLTExLTE1LCBhbmQgdGhlbiBpbmNyZWFzZXMgYWdhaW4gc2luY2UgMjAyMS0wNSwgYnV0IG5vdCBhcyByYXBpZCBhcyBpbiAyMDIwLiBDb21iaW5lZCB3aXRoIHRoZSBtb250aCBhbmQgc2hvb3RpbmcgY2FzZXMgYW5hbHlzaXMgYmVmb3JlLCBpbmNyZWFzZSBpbiBzaG9vdGluZyBjYXNlIGR1cmluZyBNYXkgaXMgY29tbW9uLCBidXQgZXNwZWNpYWxseSBoaWdoIGdyb3d0aCByYXRlIGR1cmluZyB0aGUgYmVnaW5pbmcgb2YgZXBpZGVtaWMgaXMgdW51c3VhbC4gV2hhdCdzIG1vcmUsIGl0IGlzIG5vdGV3b3J0aHkgdGhhdCB0aGUgc3ByZWFkIG9mIERlbHRhIHZpcnVzIGFsc28gYmVnaW5zIGluIHRoZSBlbmQgb2YgQXByaWwsIDIwMjEuCgpgYGB7ciBjdW1fcmVsYXRpb259CnNob290aW5nX2NvdmlkICU+JSAKICBtdXRhdGUodGV4dF9sYWJlbCA9IHN0cl9jKCJkYXRlOiAiLCBkYXRlKSkgJT4lIAogIHBsb3RfbHkoeCA9IH5jb3ZpZF9jdW1fcGVyMTAwMCwgeSA9IH5zaG9vdGluZ19jdW0sdHlwZSA9ICdzY2F0dGVyJywgbW9kZSA9ICdsaW5lcycsIHRleHQgPSB+dGV4dF9sYWJlbCwgYWxwaGEgPSAwLjgsIGNvbG9ycyA9ICJ2aXJpZGlzIikgJT4lIAogIGxheW91dCgKICAgIHRpdGxlID0gIkxpbmVhciBSZWxhdGlvbnNoaXAgYmV0d2VlbiBBY2N1bXVsYXRpb25zICIsCiAgICB4YXhpcyA9IGxpc3QodGl0bGUgPSAiQ3VtdWxhdGl2ZSBjb3ZpZCBjYXNlcy8xMDAwIiksCiAgICB5YXhpcyA9IGxpc3QodGl0bGUgPSAiQ3VtdWxhdGl2ZSBzaG9vdGluZyBjYXNlcyIpKQpgYGAKCg==